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Abstract 

The inverse /3-decay reaction, UeP — e~^n, for low-energy anti-neutrinos 
coming from nuclear reactors is of great current interest in connection with high- 
precision measurements of the neutrino mixing angle ^13. We have previously 
derived analytic expressions, up to next-to-leading order in heavy-baryon chiral 
perturbation theory, for the radiative corrections (RCs) and the nucleon-recoil 
corrections both for this reaction and for the related neutron /3-decay process. 
We investigate here the numerical consequences of these analytic expressions. 
We show that the recoil corrections are small for neutron /3-decay, but for 
inverse /3-decay, the recoil corrections are comparable in size to the RCs for 
typical energies of reactor anti-neutrinos, and they have opposite signs. It turns 
out that the RCs and the recoil corrections exhibit very different dependences 
on the neutrino energy. 
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Very recently, several experimental collaborations reported nonzero values of the 
neutrino mixing parameter, ^13 [H |21 [31 S] . Low-energy anti- neutrinos from nuclear 
reactors are well suited to determining ^13. The Double-Chooz p!], Daya Bay [2], and 
RENO |3] Collaborations have been carrying out high-precision measurements of 6'i3, 
by monitoring the inverse /3-decay reaction 

i^e+p—^e^ + n, (1) 

for the reactor-generated UeS. The accurate extraction of 613 from a measured 
positron yield in the reaction Eq.(IT]) requires a precise knowledge of radiative cor- 
rections (RCs) and nucleon-recoil corrections. An important issue in this connection 
is to what extent one can exploit the experimental data on neutron /3-decay 

n p + e" + z/e , (2) 

to have a better control of the RCs and recoil corrections to the inverse /3-decay 
reaction. 

In earlier works O El E], the relevant RCs were evaluated to order-a in the the- 
oretical framework developed by Sirlin and Marciano, to be referred to as the S-M 
approach; see, e.g., Refs. [HI [9]. Although the estimates based on the S-M approach 
are believed to be reliable to the level of accuracy quoted in the literature, it is not 
totally excluded that these estimates may involve some degree of model dependence. 
Meanwhile, the nucleon-recoil corrections for the inverse /3-decay process were eval- 
uated based on the l/mTv-expansion of the nucleon weak- interaction form factors, 
where rriN is the nucleon mass; see e.g., Refs. [6l [TO]. 

In Refs. [m [12], we proposed to use an effective field theory (EFT) approach to 
derive the RCs and nucleon-recoil corrections for the processes, Eq.([T|) and Eq.([2|). 
Analytic expressions for both the RCs and the recoil corrections to next-to-leading 
order (NLO) were presented for neutron /3-decay [TT], and for inverse /3-decay [T^ . 
Our analytic RC expressions are consistent with those obtained in the S-M approach 
by Fukugita and Kubota [5] , and by Sirlin [13] ; see also Refs. [6l [7] for earlier works. 
As for the recoil corrections, our analytic results for inverse /3-decay presented in 
Ref. [12] were found to be consistent with Ref. [TU]. It was noted, however, that 
the recoil corrections (in particular, the one arising from the weak-magnetism term) 
have very different analytic expressions for neutron /3-decay and the inverse /3-decay 
reaction. 

In Refs. [m [12] , we did not go into detailed examinations of the numerical con- 
sequences of our analytic results. The purpose of the present note is to supply such 
examinations, which we hope will shed further light on the practical significance of 
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the RCs and the recoil corrections. In this paper, after a brief explanation of the 
theoretical framework used, we present numerical results for the RCs and recoil cor- 
rections calculated up to NLO for representative observables (see below) for neutron 
/3-decay and inverse /9-decay. 

Let Q represent the typical four-momentum involved in neuron /3-decay, or inverse 
/3-decay. For neutron /9-decay, because of its low Q-value {Qp = 0.782 MeV), Q is very 
small compared with the nucleon mass; in fact, it is very small even compared with the 
pion mass. The same is true with inverse /3-decay (with Qmv-p = mp—mn—me — —l.S 
MeV), so long as we only consider cases where incident anti- neutrinos are those coming 
from nuclear reactors. Therefore, the two processes under consideration are naturally 
described by effective field theory (EFT). In particular, the well-established chiral 
perturbation theory (xPT), see e.g., Refs. provides an ideal framework 

to carry out a systematic calculation of the electroweak RCs in a model-independent 
and gauge-invariant manner. xPT treats all strong-interaction and electroweak cor- 
rections within a unified power-counting scheme. This EFT framework admits a 
perturbative expansion of relevant Feynman amplitudes in terms of two expansion 
parameters: (i) the xPT expansion parameter, Q/A^, where ~ 4:7!-/,, ^ 1 GeV 
{fir = pion-decay constant) is the chiral scale; (ii) the usual QED expansion param- 
eter, a/{2n), which governs the contributions from the electroweak RCs. We adopt 
here heavy-baryon chiral perturbartion theory (HBxPT); see e.g., Ref. [H]. In this 
scheme an additional expansion parameter, Q/m^, governs the contributions from 
nucleon- recoil corrections. Since ~ ttitv ~ 1 GeV, the Q/A^ and Q/rriM can be 
considered as one expansion parameter. Thus, the m~^-expansion is a natural part 
of our counting scheme, which dictates how to systematically incorporate the recoil 
corrections. In EFT the short-distance physics, which is not probed in low-energy 
processes, is subsumed into so-called low-energy constants (LECs). Although in prin- 
ciple the LECs can be evaluated from the fundamental theory, i.e., QCD, but a more 
pragmatic approach is to determine the LECs by fitting experimental observables. 
Once the LECs are known, we can make predictions for other measurable quantities. 
Our concern here is to carry out a HBxPT-based calculation of the RCs and the recoil 
corrections up to NLO, viz., first order in a/27r, Q/A^, and Q/niN. 

The effective lagrangian relevant to our calculation includes the relativistic lep- 
tonic weak interaction current and the heavy-baryon chiral lagrangian. To NLO, it 
is given by 

^eff = ^QED + ^NN + £^NNipi} , (3) 
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dfj^A^ — d^A^ is the elec- 



d^ + ieAfj, is the covariant derivative. The Feynman 



Cqed in Eq.(jl]) is the usual QED lagrangian, where F^^, 
tromegnetic field tensor, and 
gauge is used to fix the gauge parameter, i.e., = 1- The Cnn is the part of the 
HB^PT lagrangian that includes the photon-nucleon interactions, while Cnn^p^p is 
the part that represents the low-energy LO and NLO current- current weak interac- 
tions. Eq.(|n]) also contains the explicit forms of NLO nucleon- recoil terms dictated by 
HBxPT. Here qa = 1.267 is the axial coupling constant, while f ^ is the nucleon velocity 
four- vector, and S'^ is the nucleon spin; they satisfy v-S = It is convenient to choose 
f^ = (l,0) and = (0,a/2). In the NLO part of the lagrangian, /Xy = ^J'p—fJ'n = 4.706 
is the nucleon isovector magnetic moment. The low-energy constants (LECs), ei, 62, 
Cy and Ca, incorporate the short-range radiative physics that is not probed in a 
low-energy process. The Fermi coupling constant, Gp = 1.166 x 10"^ GeV~^, is de- 
termined from muon decay, and the CKM matrix element, \Vud\ = 0.97418 ± 0.00027, 
is given by the PDG [Hj. 



We next briefiy recapitulate the EFT-based calculations of the NLO radiative 
and recoil corrections described in Ref.[TT] and Ref.[T2] for neutron /3-decay and 
inverse /3-decay, respectively. We shall focus on experimental situations where none 
of the particle spins are monitored by the detectors. Furthermore, for the sake of 
definiteness, we concentrate on the differential decay rate for the neutron /3-decay 
process, n{pn) p{Pp) +Ue{Pu), and the differential cross section for the 

inverse /3-decay reation, i^dPu) +p{Pp) Here the four-momenta of 

particles are indicated in the parentheses; for the four-momenta, p and p^, we shall 
also use quantities defined hj p = {E,p) and Piy = {E^,Pj^). 
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For neutron /3-decay, from Ref. [TT], we have the differential decay rate 

2 



dr 



dE d(cos^e 



g{E) {l + 3g\)}C,{f3) + {l-g\)}CMf3cose, , (7) 



where /Cj(/3)'s (i = 1, 2) contain the radiative and recoil corrections, to be described 
below. Meanwhile, according to Ref. [12], the differential cross-section for inverse 
/3-decay is given by 



da., 



inv—fj 



GpVud 



f{E) {l+3gl)G,{P) + {l-gl)G^{P)Pcose, 



d(cos^e) V V2 

where ^j(/3)'s {i = 1,2) contain the radiative and recoil corrections, see below. In 



Eqs.(I7j) and (IHl), (3 =\p\/ E =^ E'^ — ml/ E is the velocity of the outgoing electron/positron 
and cos(6'e) =Pu P- The phase-space factor, g{E), that includes the (9(m^^) correc- 
tions is given by [H] 



9iE) 



F{Z,E) E'^l3{E'^'''-'-Ef 
2^^3 



1 + (3^ 3/3Ecos^e) + 0{m 



rriN 



N 



, (9) 



where E"^"'^ = (m^ — -|- Tn1)/2mn = E^ + E + 0{m]^^) is the maximum or "end- 
point" energy of the emitted electron in neutron /3-decay, and F{Z,E) ^ l + air/ (5 
{Z=l) is the usual Fermi-Coulomb function. The phase-space factor, f{E), in Eq.® 
was obtained in Ref. [12] a^ 



fiE) 



TT 



1-— fe-^cosd+0(m^') 



(10) 



The above velocity-dependent functions, /Ci(/3) and ^i(/3) (z = 1, 2), contain terms 
coming from the RCs and those from nucleon-recoil corrections. 



Q,[fi) = i + -^e;™'^(/3) + ^e;[--'(/3). 



(11) 

(12) 



Note that for inverse /3-decay, the positron energy, E, and velocity, /3, contain terms of ©(m^^), 



e.g., = E{1 — 0{mj^ )), where E = E^, — (m„ — mp); see the discussion leading to Eq. (14) in 
Ref. ini. 
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As for the recoil corrections, we note that, since the kinematical l/mj\f corrections 
are already included in the phase-space factors, g{E) and f{E), the recoil corrections 
in Eqs. fllll) . f|T2l) are dynamical ones coming from the 0{m]^^) terms in Eq.([6]). 

The HBxPT calculation, up to NLO, of the RCs was given in Ref. fllj for neutron 
/9-decay, and in Ref. [12] for inverse /9-decay. Here we simply quote the results, 
referring the reader to Refs. [TT] and [12] for the details of their derivationjf] The 
results are given as: 
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(13) 
(14) 
(15) 
(16) 



The so-called "inner" RCs, which are independent of /3, are encoded in the LEG, 
e^(/i2). The "outer" RCs, 6'i\(3), 5i^H/3), ^out{l^) and 5o«t(/3), constitute the well- 
known, model- independent, long-distance QED corrections that do not contain any 
hadronic effects. Analytic expressions for these outer corrections can be found in, 
e.g., Refs. [m 112]. Note that not only the ultraviolet-divergent and scale dependent 
terms are subsumed in the LEG, ey(/i^), but also all the infrared- divergent terms of 
0{a) are simultaneously canceled, leading to finite final results. 

As for the recoil corrections, /C[^™*'(/3)'s {i = 1, 2) appearing in Eq. (fTT]) (pertaining 
to neutron /3-decay), and ^[^™*'(/3)'s {i = 1,2) appearing in Eq. (fT2|) (pertaining to 
inverse /3-decay rate) are given by (see Refs. [TT| [12]). 
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recoil ^l^'j 
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gl 



(17) 



^ It is to be remarked that the HB^PT rcsuhs are consistent with those obtained in the the S-M 
approach [H E [Zl |8] . 
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A noteworthy point in the above expressions is that the terms involving 2gA^y ~ 
12 arising from the "weak-magnetism" interactions in Eq.([H]) constitute the most 
dominant recoil corrections. 

The above formulae allow us to take a formally consistent account (up to NLO in 
HBxPT) of the RCs and recoil corrections for both neutron /9-decay and the inverse /3 
process. Since these two processes involve the same LEG, e{^(/i^), as the only unknown 
parameter, it is in principle possible to use neutron /3-decay data to fix ey(/i^), and 
make model-independent predictions for the inverse /3-decay cross sections. This 
statement is true to the extent that the other quantities which go into the above 
formulae and which are usually treated as "known" quantities are indeed known with 
high enough precision. To what degree the existing uncertainties in the values of 
gA and the neutron lifetime affect the outcome of such an analysis is a subject that 
deserves a careful study. Relegating this investigation to a future study, we choose 
here to use an estimate for ey(/i^) deduced in Ref. pjj by comparing the results for 
neutron /3-decay obtained in HB^PT with those obtained in the S-M approach [H [9] . 
This comparison yields 

~eR[^^=nil)^A\J'^+\J'^+2C+A^ = 19.5 ±0.7, (19) 

where the following input has been used. From a "liberal" choice of the Ai-resonance 
mass, one has myi~1.2±0.6 MeV; the electroweak "loop" -function, C, is estimated to 
be C = 0.829(l±0.1), while a perturbative QCD calculation gives Ag = —0.34 (see, e.g., 
Ref. [IH])- Thus the size of the RC involving the LEG in Eqs. (fT3l) ^(IT6l) is estimated 
to be (a/27r) ey(m^) ~ 0.023. This short-distance (inner) contribution to ey(m^) 
is dominated by the well-known box diagrams involving Z-boson exchanges [HI |9]. 
The numerical results reported in what follows have been obtained with the use of 
e^{lj.^=ml) =19.5 ±0.7, 

To illustrate interplay between the RCs and recoil corrections, and to exemplify 
the differences between the recoil corrections for /3-decay and inverse /3-decay, we con- 
centrate on the angle-integrated observables, viz., the differential decay rate, dTis/dE, 
for /3-decay, and the total cross section, amv-p^ for the inverse /3-decay reaction. This 
means that we are concerned with /Ci and Qi in Eqs. ([7]), (jS]), or equivalently, with 

^rad^ grad^ J^recoil ^ grecoil EqS.([IID, dH- 

In Fig{T] we plot ^/C™*^ for neutron /3-decay as a function of the anti-neutrino 
energy, E^, {0 < < 0.78 MeV). The pair of RC curves in the figure indicates an 
error-band associated with the above-mentioned uncertainty in ey(/i^ = m^). Figj2] 
shows ^Gl"'^ for the inverse /3-decay reaction for a range of E,^ corresponding to 
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reactor-generated anti- neutrinos {E,^ < 8 MeV). Here again, the pair of RC curves 
represents an error-band reflecting the uncertainties in ey(/x^=m^). 

FigJUalso shows "^r'^"*' /3-decay, while Figl2] includes for the in- 

verse /9-decay reaction. It is seen that /C?['^'^°*' is much smaller than ^[^™*'^ which is 
not unexpected, because E,y for reactor neutrinos is significantly higher than the anti- 
neutrino energy involved in neutron /3-decay {Ei,< 0.78 MeV). We point out, however, 
that the difference between /C^*^™*' and ^j^'^™*' is further enhanced by the fact that, 
as seen in Eqs. fll7l) and f[T5]) . there is near cancellation among the various terms 
contributing to /C^^™*', whereas they add up for 

It may be of interest to consider the "total recoil" correction, which represents the 
combined effect of the m]^^ terms in the lagrangian, Eq.(|6]), and the correction 
in the phase space factor, Eq.Q. The differential decay rate, Eq.Q, after angle- 
integration can be rewritten as 



dE ~ dE r + 27r^i ^'^ mj. 



1 + — /C™'^(/3) + /cr"'-'"'='=°''(/3)l , (20) 

2n rriAT J 



where dT^ /dE = {GFVud)^[E{Z = 1, E)E'^f3{E""''=-EY /{2n^)]{l + 3 g\) is the stan- 
dard LO result, and the total recoil correction, /(;;t°t'^i-''*^coii given by 



^total-recoil = /C;'"^°*' (/3) + 3E - E™"^ - ( ^^T^ 1 ^ + ^i^N^) ' (21) 



Meanwhile, Eq.(|8]) (after integrated over the angle) can be rewritten as 

cr.n.-p = aJ°L, {l + £^r'(/3) + 1(/3)} , (22) 

where (T-°^_^ = (G'irK(cz)^(-^^/3/vr)(l + 35'^) is the standard LO result, and the "total" 
recoil correction, gf^^^^-^^^^^^p) ^ is given bj0 

^total-recoil(^) = gi-\P)-E- (^)(^. +^^) + (tTm)^^ + ^^"'^'^ ' 

(23) 



^ The last term in Eq.(20) in Ref. [T^] involved an error in its evaluation. In Eg. ([25]) in the 
present article, this error has been corrected. 
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where E = Ey — {mn — mp) = E^ — is the "lowest-order positron energy" {i.e., the 
energy corresponding to a case where the neutron recoil is neglected), and /3 is the 

corresponding velocity: /3 =\J E"^ — m1/E% 

FigJD displays /C!°*"^-^"™^ while Fig J2] shows ^total-rccoil^ ^^^^ ^^^^ |^total-recoil| 

is significantly smaller than j^totai-recoiij^ comparison between the RCs and the 

recoil corrections, FigiHand Figj2] indicate that the recoil corrections are much smaller 
than the RCs for neutron /3-decay. By contrast, for inverse /3-decay process, Q\°''^ and 
^total-recoil Comparable size and have opposite signs. It is also to be noted 

that, as seen in Fig. [21 for the inverse /3-decay reaction the recoil corrections and the 
RCs have very distinct energy dependencies. Thus both corrections must be carefully 
considered in analyzing high-precision data used to deduce ^13. 

To summarize, we have presented a brief discussion of numerical consequences 
of the radiative corrections (RCs) and the recoil corrections for neutron /5-decay, 
Eq.(l2]), and inverse /3-decay, Eq.([I]), calculated up to next-to-leading order (NLO) in 
HB^PT. The numerical results reported here are obtained with the use of the value 
of the LEC, ey(/i^), that has been deduced from comparison of our HB^PT results 
with those obtained in the S-M approach. This hybrid nature of our analysis should 
eventually be replaced by a more rigorous determination of ey(/i^) through a direct 
fit to the neutron /3-decay data. It is expected, however, that the basic features of 
the RCs are already visible in the present hybrid treatment. The recoil corrections, 
which up to NLO do not involve the LEC, are not affected by the hybrid nature of 
our numerical calculation. 
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Figure 1: For neutron /3-decay, the radiative correction, ^ /C™'' in Eq. (fT3|) . and 
the recoil correction, /C^^™*' in Eq. ( !T7|) . are plotted in % as functions of 
the emitted anti- neutrino energy Ey. The dashed curve labeled "Total Recoil" gives 
_j_ ^total-recoil Eq. (|2Tl) . The pair of radiative correction curves gives an "error 

band" that reflects uncertainties in the values of the LEG, ey(m^), see the text. 
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Figure 2: For the the reaction z/g + p — )• e"^ + n, the radiative correction, ^ in 
Eq. ( !T3|) . and the recoil correction, m^^ grecoii -g^ (|T8|) . are plotted (in per cent) as 
functions of the incoming anti- neutrino energy E,^. The dashed curve labeled "Total 
Recoil" gives _i_^totai-rccoii Eq.(l23ll. The pair of radiative correction curves gives 

an "error band" that reflects uncertainties in the values of the LEG, ey(m^), see the 
text. The vertical dotted line indicates the threshold {E^c:^1.8 MeV). 



11 



